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ABSTRACT 

We develop a perturbative expansion which allows the construction of non- 
abelian self-dual SU(2) Yang-Mills field configurations on the four-dimensional 
torus with topological charge 1/2. The expansion is performed around the 
constant field strength abelian solutions found by 't Hooft. Next to leading 
order calculations are compared with numerical results obtained with lattice 
gauge theory techniques. 
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1 Introduction 



Self-dual Yang-Mills fields are fascinating mathematical objects that play an 
important role in both Physics and Mathematics. In the Physics literature 
they emerged through the introduction of the BPST instanton [|l|, the mini- 
mum action configuration in the sector of topological charge Q equal to one. 
This triggered the joint effort of physicists and mathematicians in the search 
for multiinstanton configurations. This had its reward with the ADHM for- 
malism p[, a general set up for the construction of self-dual configurations 
with vanishing field strength at infinity (which, from the topological point of 
view, corresponds to fields living on a four-sphere). However, there are cer- 
tain instances in which other types of boundary conditions are relevant. For 
example, in considering finite temperature field theory, one is interested in 
studying configurations which are periodic in thermal time. This periodicity 
might also be used as a device to study monopole-like objects. In the same 
fashion, additional periodicities might have other uses and interpretations, 
as argued in Ref. 0. This justifies an initial interest in the study of self- 
dual configurations which are periodic in all euclidean space-time directions. 
Geometrically, it corresponds to the study of self-dual gauge fields on the 
torus. It might be considered as the next step after the case of gauge fields 
on the sphere. The study of gauge fields on the torus (see |Q for a review of 
this topic) brings in a new topological richness, whose appearance, physical 
interpretation and usefulness was put forward by 't Hooft p. In addition, 
periodic self-dual configurations constitute a simple mathematical example 
of dense multiinstanton configurations, which, as advocated in Ref. 0, might 
turn out to be a better description of the confining vacuum than other dilute 
multiinstanton pictures. 
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The construction of self-dual configurations on the torus has met very 
limited success. For particular values of the torus sizes, there is a class of 
solutions which is known since early times 0. These configurations have 
constant field strength, and they are, in some sense, of an abelian nature: 
all the spatial components of the electric and magnetic fields are parallel 
in colour space. For other torus sizes, we know, through numerical tech- 
niques P, H, 0, that the solutions are quite different. They turn out to be 
lumpy and very non-abelian: different electric field components are mutually 
orthogonal in colour space at the center of the lump. Fortunately, there has 
recently been substantial progress in studying the mixed situation of instan- 
ton configurations which are periodic in only a few directions. In particular, 
in the case of only one compactified direction, corresponding to the physi- 
cal situation of finite temperature, the most general solution of topological 
charge one, the caloron, has been found [0. Considerable progress has also 
been achieved for the case of gauge fields on x i?^ (doubly periodic in- 
stantons), by the work of mathematicians ||ll|] and physicists [0, |T2|, [T^. In 
these developments a crucial role is played by the Nahm transform, a duality 
transformation which maps self-dual configurations on tori with dual sizes 



14, 15 . For the T x R case, the calculation of the abelian Nahm dual of the 



charge one instanton [jT6| represents a step forward. The caloron, the doubly 
periodic instantons and the x R instanton can be viewed as limiting cases 
of configurations on the torus where either three, two or one directions are 
taken to be very large with respect to the others. 

This paper is a step towards an analytical understanding of non-constant 
field strength self-dual configurations on the torus. The strategy is to consider 
torus sizes which depart only slightly from those in which there exists a 
self-dual constant field strength solution. Our construction is based upon 
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perturbing around the constant solution (at a torus size in which it is not self- 
dual), and imposing self-duality to the resulting configuration. A systematic 
perturbative expansion arises that allows the construction of the self-dual 
solution. We show that the solution exists order by order and is unique 
up to gauge transformations and space-time translations. This is done in 



section 2. Our approach is intimately related to the study of van Baal |jT7 
who considered perturbations around constant field strength solutions. 

In section 3 we proceed to compare the lowest non-trivial order results 
obtained from our expansion with the exact solution obtained by numerical 
methods. This serves to quantify the rate of convergence of the series for 
various torus sizes. Finally, in section 4, we investigate the interplay of our 
perturbative expansion with the Nahm transform. The paper is closed by 
section 5, where the conclusions and possible extensions are presented. 



2 The construction 

Let us consider SU{2) gauge fields living on a torus of size Iq x li x I2 x l^. 
Under a translation by one of the periods, the gauge potentials and fields 
transform by a gauge transformation: 

A,{x + e^) = [Q^{x)]A,{x) , (1) 

where is a 4-vector of length along the fi-th direction, and are the 
twist matrices. The compatibility conditions of the previous equations are: 

^f^ix + e^) VL^{x) = exp{mn^^} Q^{x + e^) Q^{x) , (2) 

where the elements of the antisymmetric twist tensor 12^1, are integers defined 
modulo 2. In what follows we will choose nos = "^-12 = —"".30 = —n2i = 1, 
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with the remaining components being zero. For the twist matrices we will 
take 

Q^{x) = exp{?^n^i.y^r3} , (3) 

which is consistent with our choice of twist tensor. The symbols label the 
Pauli matrices. 

For torus sizes such that loh — hh, there exist self-dual configura- 
tions satisfying the previous boundary conditions and having constant field 
strength. What we will do is to consider a slight deviation from this situation 
controlled by the parameter: 

yV 



where V = klikh is the torus volume. Without loss of generality wc can 
assume that A is positive. In this case there exists a constant field strength 
configuration with vector potential: 

Bn{x) = —n^u^Ts . (5) 

This gauge potential gives rise to a field strength of the form G^^t^, where: 

G,. = 7rff , (6) 

The only non-zero components are Cos and G12, which become of equal 
magnitude at A = 0, rendering the solution self-dual. With our choice of 
twist the constant field strength configuration has topological charge Q — 1/2 
and, for A = 0, total action 47r^. 

Now, let us consider perturbing around this gauge potential: 

A^{x)^B^{x) + S^{x)t^ + W^{x)t+ + W;{x)t_ , (7) 

where wc have decomposed the additional field into different colour compo- 
nents. The matrices t± — |(ri±zT2) are standard. The boundary conditions 



on the gauge fields translate into the real function S^{x) being periodic on 
the box, and the complex function W^{x) satisfying 

Wp{x + e^) = exp{OT rif,^^} Wp{x) . (8) 

We will make the following gauge choice, consistent with these boundary 
conditions (the background field gauge): 

dpAp{x)-i[Bp{x),Ap{x)]=0 . (9) 

We will now demand that the resulting gauge field is self-dual 

F^,(x) -F^,(x) = 0, (10) 
F/^f (3^) = 2^^ll'paFpa{x) , with eoi23 = 1 , 

which will be interpreted as equations for the functions Sp{x) and VF^(x). 
The best way to express these equations, together with the gauge fixing 
condition, is to use the matrices cr^ = (I, —tr) and = (I, tr) = crjj. These 
matrices satisfy: 

^^iCr^ = fj'^^aa , (11) 

where r]'^^ is the 't Hooft symbol, such that r/°^, = 5^1, and the r/^^ are a basis 
of the antiself-dual tensors. Now contracting Fp^{x) with a^a^, we project 
out the self-dual part. Hence, we might rewrite equation ([ToD as follows: 

OS = -|G + ^{W^Wc - W^W) (12) 
DW = t{S^W - W^S) , (13) 

where S = Sp{x)a^ and W = Wp,{x)ap are 2x2 matrices, S\ their 
adjoints and the parameter A is equal to 1. The matrix Wc is the charge 
conjugate: 

Wc = T2W*T2 . (14) 
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The matrix G = Gni,au(Ju is given by: 



G = 2m^T^ . (15) 



It vanishes when A = 0, since then the constant field strength configuration 
is self-duaL Finally, we define the following differential operators: 

(16) 
(17) 
(18) 
(19) 
(20) 



d = 




d = 




D = 




D = 




D, = 


df, + m 



Since at A = the correction terms W and S vanish, we can think of 
A as a perturbation parameter. Rigorously speaking this is not quite so, 
because A depends on the torus sizes, and they enter also in the boundary 
conditions. To keep a truly perturbative parameter in our expansion we 
introduced in Eq. ([T^) the parameter A, whose interpretation will be clear 
later. Our goal will be to solve Eqs. (|T^-(|T3p for arbitrary values of A, as 



a perturbative expansion in powers of AA. For that we have to expand the 
unknown matrices S,W in powers of V AA. We see that the equations are 
consistent with W carrying odd powers and S even powers: 

W = -^^(AA)'=W^('=) (21) 



k=l 

oo 

S = ^(AA)'=5('') . (22) 

k=l 

In the following paragraphs we will show that it is possible to solve the set 
of equations (|T^-(|I3|) order by order in V AA. Finally, setting A = 1 one 



recovers the solution of the self-duahty condition. However, even for this 
value higher orders in the expansion are suppressed by powers of Thus, 
as we will verify later, we expect the first few terms of the expansion to 
approximate the self-dual solutions for small values of A. 

On the other hand, the solution for arbitrary A can be interpreted as the 
solution of the following modified self-duality equation: 

F^,{x)-F,,{x) = {l-X){G^,-G^,)r, . (23) 

For A = the constant field strength configuration is a solution, while for 
A = 1 we recover the self-duality equations. 

Now let us address solving the equations order by order in A. Notice first 
that W^-^^ satisfies the equation: 

W(^) = . (24) 

This equation has non-zero regular solutions, as can be deduced from the 
index theorem. As we will see later the general solution has the form: 

= ^(x) , (25 

\ j 

where K^^^ and Q^^^ are two arbitrary complex numbers and \I'(x) is a func- 
tion whose explicit form will be given later. Having seen that Eq. (|13D has a 
solution ioT k = 1, let us address the question of whether it also has a solu- 
tion for all values of k. For that purpose one has to investigate the adjoint 
of D, and see whether its kernel is null or not. Indeed, it is easy to see that 
the kernel of the adjoint vanishes, and hence Eq. (0) has a regular solution 
no matter what the left-hand side is, provided it is regular. Explicitly, one 
such solution is given by: 

M^W = L)f/W (26) 



k-1 



where DD f/C^) = z ^ [S^'^-^^^W^^^ - W^i'^t^^'^-')) . (27) 



1=1 



The existence of a unique solution U^'^^ to Eq. (^7]) follows from the invert- 



ibility of the operator DD. In terms of this particular solution, and using 
(p^)-(p5|), the most general solution of Eq. (|T3|) is given by: 



where K^''^ and Q^''^ are complex constants. 

Now let us study the solution of equation ( |T2l) order by order in A. It is 
easy to see that both the left as the right hand sides are periodic functions in 
the box. They can hence be expanded in Fourier series. The solution can be 
obtained by equating the corresponding Fourier coefficients. However, notice 
that the left-hand side has no constant term. Hence, if the right-hand has a 
non-zero constant Fourier term, the equation has no solution. This can be 
expressed more formally by saying that the kernel of d is nontrivial. What 
we will now show is that in solving the equation for S^''^ it is possible to 
fix the constants K^''^ and Q^''^ appearing in the solution to (|l^) to order k, 
by the condition that the lowest Fourier component of the right-hand side 
vanishes. For that we have to explore the effect of the replacement ( pH]) in 
Eq. (|1^) to order (A A)'''. The dependence of the right-hand side of Eq. ([12D 
on K^''^ and Q*-*'-' to this order, is contained in the following term: 



where Cg = 23ft(Q(^)gW* - i^Wi^W*) and c+ = 2(irW*QW + ;r«*Q('=)). 



The symbol 3? denotes the real part of its complex argument. By choosing 
]^(k)^ Q(k) appropriately the term within parenthesis can be made equal to an 




(28) 



^(x)l^ (03X3 + c+r+ + c;r_) 



(29) 
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arbitrary hermitian, traceless 2x2 matrix. Since the lowest Fourier coefficient 
(the constant one) of is non-zero, the constants can be chosen such 

that the whole right hand side of Eq. ([I^ ) has vanishing constant Fourier 



term. Actually, this fixes 3 of the 4 real parameters which enter K^''\ Q^^\ 
The remaining one corresponds to the symmetry associated to global colour 
rotations in the 1 — 2 plane, which leave invariant (we will comment upon 
this property on the next paragraph). 

Having shown that the solution of our set of equations exists order by 
order in our expansion in V AA, we have now to analyse uniqueness. Indeed, 
on general grounds we know that the solution is non-unique. This fact is 
associated to the existence of transformations which change one solution into 
other. We already mentioned one: global gauge transformations of a certain 
kind. These are the residual gauge transformations that are not gauge fixed 
by Eq. (^. They are associated with the freedom to multiply any solution 
matrix W hy a. constant phase. The other transformations are space-time 
translations (followed by an appropriate gauge transformation to preserve 
the gauge fixing condition). This latter symmetry manifests itself under the 
form of a non-uniqueness for the solutions of Eq. (|T^) to any order k: notice 
that we are free to add an arbitrary constant matrix to S in the left-hand 
side of the equation, which would entail the four real parameters associated 
to a translation. 

The best strategy in solving the equations is to fix a unique solution by 
constraining these transformations. This we will do by imposing the following 
additional conditions: 

^{Wi2{x = Q)) = \Wi2{x = Q)\ (30) 
jdxS{x) = Q , (31) 
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(If Wuix = 0) = we take ^{Wn{x = 0)) = \Wu{x = 0)|). It is now 
completely clear that the procedure leads to a unique solution order by order 
in V^. 

It is useful to derive expressions for the field strength tensor itself. Just 
as we did for the vector potential we might expand in colour components: 

= Fj^fix) r3 + F(+)(a;) r+ + F;y\x) r_ . (32) 

Now since the field is self-dual, we might contract it with the matrices a^jWy 
to obtain a traceless hermitian matrix combining the three spatial directions: 

^^^^ = ^F(t\x)a,a^ 

= Y^W^'l^^^c-WS^) (33) 



^^'^ - -{Fl:t\x)a,a^ 

= l(^-^^^)rs-'fS^-\{WW^-W.W^J . (34) 

Having set up the full procedure for calculating the potentials and fields 
in powers of A, let us now exemplify it by computing the first terms in this 
expansion. These results will be used in the next section. The starting point 
is the equation for W^^^ (Eq. (p^)). We mentioned previously what is the 
form of the solution. Let us for the moment skip the proof and also the 
determination of "^{x) and proceed. The next step is to look at the equation 
for 5*'-^''. As mentioned in the general case, both the left and right-hand sides 
can be expanded in Fourier coefficients. The condition that the constant 
coefficient of the right hand side vanishes imposes a constraint on K^^^ and 

Q«*ir«=0 (35) 
|g(i)|2_|^(i)|2^^ (36) 
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where we have fixed the normahsation of "^{x), such that its constant Fourier 
coefficient is equal to one. The previous equations lead to 

ir« = 0; Q(^) = ^ , (37) 



where we have used (pOf ). 

Now the equation for S^^^ reads: 



dS^'^ = ^m^)\' - 1)t, . (38) 



This can be solved together with Eq. (pID to give: 



^W = ^(9%3 , (39) 



where h{x) is a periodic function on the box, solution of the equation: 

□/i(x) = |^(x)|2 - 1 (40) 

and □ is the 4-dimensional Laplacian. The previous equation can be solved 
by expanding both sides in Fourier series and equating. 

Let us now work out the details of the solution to Eq. (j2^). For future 
purposes we will consider a more general equation: 

D~gip{x) = , (41) 

where is a two component vector. The operator Dg is given by: 

D-, = a,id, + t7rq^) , (42) 

where g is a constant. The vector of functions (p{x) is required to satisfy the 
boundary condition: 

X 

ip{x + e^) = exp{t7T q n^^-^} <^{x) (43) 

''1/ 
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which is only consistent for integer q. It is easy to see that the operator Dq 
preserves this boundary condition. 

Now it is seen that locally a solution of Eq. (^) takes the form: 

(^g(a;) k+(mo,Ui) 
where we have introduced complex coordinates: 



(44) 



u^ = ^{x^ + i Uf^^Xy) (45) 

and u* are the complex conjugates. These coordinates are not independent, 
and satisfy: 

M/i = ■J-IT-ij.uUju ■ (46) 

We might for future benefit introduce the complex constants: 
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The function <fq is given by: 

ipq{x)=exp{-^{xl + xl)-^{xl + xl)} . (48) 

The boundary conditions Eq. (^) impose constraints on the value on the 
holomorphic and anti-holomorphic functions k±: 

{u + -) 

K+(x + e„) = expjvrgz — —} k+(x) (49) 

{u* + -) 

K_(x + e„) = exp{— 7rg« — —} k_{x) . (50) 

Choosing uq and ui as our two independent complex variables, we might 
write: 

k+(mo,Mi) = exp{7rg2(— 2- + —^)} k+(mo,Mi) (51) 

Ztq Zti 

/«_(«*,«*) =exp{-7rgz(-^ + -^)} K_(n*,<) (52) 
13 



The functions k± are periodic in their arguments with period 1 and satisfy: 

K,+ {uo + To, Ml) = exp{-27rgzMo - '^q^To} «5+(mo7 ^i) (53) 
K+{uo, Ui + Ti) = exp{— 27rgmi — nqiTi} k^{uq, Ui) . (54) 

For q = 1 these are precisely the conditions satisfied by the Riemann 6 
function |[TB|. Actually, up to a multiplicative constant, this function is the 
only (regular) holomorphic function satisfying these boundary conditions. 
Similarly, one obtains that the equation for k_ has no regular solutions. 
Hence, for the q = 1 case we have arrived at the solution given in Eq. (p5D, 
and determined the expression for the function \E'(x): 

^(^) = \ ~rr 6^p{-7-r(^3 - ^^^s^^o) - rrixj- tx2Xi)} e{uo, tq) 6{ui, n) 

V 'on 'o'3 nh 

(55) 

The multiplicative factor preceding the right hand side of the previous ex- 
pression is determined by the condition that the lowest Fourier coefficient of 
|^(x)p is unity. 



3 Comparison with numerical results 

In the previous section we have set up a general procedure to construct the 
form of the gauge potentials and field strengths for SU{2) self-dual solutions 
on the torus with twist tensor rioa = "^12 = 1- The result is an expansion in 
powers of \/A, where A is defined in Eq. (^. However, we do not have an 
analytical estimate of the size of the coefficient. Our purpose in this section 
is to test the rate of convergence of the expansion by comparing the results 
obtained from the first non-trivial order with the exact result as obtained by 
numerical methods on the lattice. 
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We will restrict to the analysis of the gauge invariant traces Tr(EjEj), 
where Ej = Fo, are the electric fields. To lowest order in our expansion, 
we fall into the constant field strength configuration, and the only non-zero 
gauge invariant trace is Tr(E3). The next correction is order A and vanishes 
for Tr(EiE2). It also predicts that Tr(E2) = Tr(E^). Making use of the 
general formulas (|33|), (0) and substituting the explicit form of S^^^ and 
yj/(o) (^Eqs. (^9]), (0)) one arrives at the following result, valid to order A: 



Tr(E^(a;)) = Tr(E^(x)) = A-= \Do^{x)\' , (56) 



Tr(Ei(x)E3(x)) = -An' (^-^^y^) + ^'^'^ ^^^^ ' ^^^^ 

2 / ^0^3 + hh 



Tr(E2(x)E3(x)) = Avr^ [ ys/2 ) i^odi - d,d,) h{x) , (59) 
where and h have been defined in Section 2 (Eqs. ( [55| ) and (P0|)). Using 



the standard representation of Riemann's theta function [18 



9{u, t) = exp{27r«nM + vrm^r} (60) 
nez 

one can easily obtain the Fourier coefficients of all the functions appearing in 
Eqs. ( p^ - |59| ). As for the numerical comparison, summing the first few hun- 
dred terms of the Fourier expansion allows to compute these functions with 
negligible errors. It is also extremely simple to use these Fourier coefficients 
to integrate analytically over some of the four real coordinates, to arrive at 
a quantity better suited for graphically displaying the comparison. 

A numerical approximation to the exact solutions of the self-duality equa- 
tions, with which the results coming from the perturbative approximation 
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that we have developed are to be compared, can be constructed by means 
of standard lattice gauge theory techniques ^. We will use for this pur- 
pose an e = overimproved cooling procedure, that was found in previous 
works 1^ to be able to produce very accurate approximants to continuum 
self-dual fields. In particular, it allows to extract the exact values of the 
gauge invariant densities Tr(EjEj) under concern up to 0{a'^) corrections, a 
being the lattice spacing (whose precise definition we will discuss below). 

To explore the accuracy of the next to leading term in the perturbative 
expansion with varying values of the remaining parameters, we will consider 
tori of lengths (/q = + = h,h = h,h = h)- The results will then 
depend on the perturbative parameter A = e/y/1 + e, and on the ratio Is/k 
measuring the degree of spatial asymmetry of the torus. Thus, by keeping 
< e ^ 1 we will remain within the perturbative regime. On the other 
hand, we will vary the ratio Ig/h between 1 (the more symmetrical case) 
and 0, were the 6 functions and their derivatives entering into the analytical 
expressions are well described by polynomials times gaussians. 

In each numerical solution is obtained in a lattice with a number of 

points = l^/a along direction fi. To define the lattice spacing a, we need 
to fix a unit. In our case we take /1/2 = 1- This is justified by noticing from 
the expression for \1/ in Eq. (^) that the region having nontrivial structure 
in the action density at nonzero A is of size y/hh- 

We will present the results of the comparison of the analytical results 
Eqs. (^^ ^U]) with their numerical counterparts, for three different configu- 
rations, having different values of A and Ig/h- The lattice sizes that we will 
use, together with their associated A and Ig/h values, are detailed in Table |l|. 

In Figures |l|, § and ^ we show the numerical and perturbative results 
for the integrated electric field densities ^f^^ {xo,xi) = J dx2da;3Tr(E|(x)) 
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Lattice 


Size 


A 




A 


13 X 12 X 12 X 12 


0.080064 


1.00 


B 


21 X 20 X 8 X 8 


0.048795 


0.40 


C 


41 X 40 X 6 X 6 


0.024693 


0.15 



Table 1: Lattices used in the comparison, and their associated A and Is/h 
values. 

and <l>fi^(xo,Xi) = $g^(a;o,Xi) = j dx2dx-{£i{Y?^{x)), and $g^(a;o,Xi) = 
j da;2dx3Tr(E2(x)E3(a;)), respectively (notice that ^^g"* vanishes at the present 
perturbative order, despite the fact that Tr(Ei(a;)E3(a;)) does not, because of 
the particular form of the expression for this latter quantity). The qualitative 
agreement is clearly good. The main features of the exact solution are present 
in the analytical expression. It is possible to obtain a graphical quantitative 
measure of the comparison by integrating the previous densities over an addi- 
tional coordinate to yield the time profiles $33^ (xq) = / da;ida;2dx3Tr(E3(x)) 
and $[^''(a;o) = *&22^(a;o) = / da;ida;2dx3Tr(E^(a;)) (similarly to what hap- 
pened with $13^ $23"* vanishes despite $23^ does not). The comparison for 
these quantities is displayed in Figure ^. 

Let us briefly comment some salient features of the solution. As for the 
density of the component E3 of the electric field, a hole appears overlaying the 
flat background supplied by the zero-order constant abelian field. The width 
of this structure is, as we mentioned above, proportional to \/Tsk, as can be 
derived from an analysis of the perturbative expression for the potential, and 
its contribution to the total action is of order A (cf. Eq. (^)). Meanwhile, 
the action density associated to the other components of the electric field, 
whose contribution to the total action is as well of order A in the perturbative 
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approach, exhibits in the asymmetric torus case (configurations B and C) 
a double lump structure, again of size \/lsh, with the maxima aligned in 
Euclidian time. We see that the perturbative approximation is most accurate 
for the configuration A, despite the relatively large value of A associated to it. 
This fact indicates that the convergence behaviour of the perturbative series 
depends on the value of the asymmetry parameter Is/h, in such a way that 
it is worse the more asymmetric the torus is chosen. The overall conclusion 
is, anyway, that for values of A in the probed range, between 0.02 and 0.09, 
the NLO perturbative result constitutes a good approximation to the exact 
solution. 

Once the convergence behaviour of the perturbative series for small values 
of A has been checked to be good, it would be also interesting to study to 
what extent the NLO result at hand remains useful to describe solutions 
occurring at larger values of A. This possibility is tempting because it would 
open the door to apply our results to improve the analytical control over 
some particularly interesting fields. For instance, it is known ^ that in a 
torus of geometry It x with Ig/k ^ 1 the solutions approach self-dual fields 
on X R, the approximation being already remarkably good for It ~ 3/^. For 
the considered twist and SU{2) gauge group a.Q = 1/2 solution is obtained, 
whose action density displays a single lump exponentially decaying in the 
large direction of length /j, and whose width is controlled by Ig- In this 
geometry, and setting It = /^(l + 5), one has A = 6/^/l + 6, and the analysis 
would proceed by moving from the case 6 = 0, where the torus is symmetric 
and the solution is the abelian one, to values of 5 ~ 1, where the features 
of the X R solution would start to arise. Having performed computations 
on lattices of sizes Lt x L^, with Lg = 12 and Lt ranging from 13 to 48, we 
have found that the self-dual configuration evolves smoothly with changing 
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S. Unfortunately, the NLO perturbative approximation begins to deviate 
substantially from the exact result before the interesting regime 5 > 1 is 
reached. 

In the same spirit one could investigate other torus geometries, e.g. x 
with Is/ It <^ 1, which in some cases is known to lead to limiting x i?^ 
solutions with a vortex-like structure 0, or k x 1^ with Ig/k ^ 1, which 
leads to the x caloron solutions [1^. In these cases, we would expect 
a similar behaviour to that found for the x R case. 



4 Nahm transform 



Nahm's transformation [14, 15] maps self-dual configurations on the torus 
onto other self-dual configurations. The modifications necessary to cope with 
twisted boundary conditions have only been worked out recently I^O], |21|] . In 



general, the transformation changes the twist tensor and torus sizes and maps 
the rank of the group (A^) and the topological charge (Q) onto each other 
through the formula: 

Q^Q' = N/No (61) 
N—^N' = QNo , (62) 

which preserves the dimensionality of the moduli spaces {QN = Q'N'). The 
integer constant Nq depends on the twist. This transformation provides an 
interesting tool for studying (anti-)self-dual gauge fields on the torus. 

In previous sections we have expressed certain self-dual potentials as an 
expansion in the parameter A. It is henceforth interesting to analyse the 
interplay of this result with the Nahm transform. First of all we should find 
out the general properties of the Nahm transform for the configurations in 
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question. In our case the group is SU{2) and the configuration has nontrivial 
twist tensor uq^, = = 1. This imphes that the parameter Nq = 4. Fur- 
thermore, the topological charge Q of these configurations is determined by 
the twist matrices and equals Q = |, as for the corresponding constant field 
strength configuration. Hence, according to the formulas given above, the 
Nahm dual is again an SU{2) solution with topological charge Q' = ^. Now 



we can make use of the results of Ref. [201 to determine the twist tensor and 



torus size of the Nahm transformed field. Indeed, the Nahm dual twist tensor 
is equal to the original one, and the torus size is given by: ^ >^ ^ >^ ^ >^ 
Thus, except for the different torus size, the Nahm transformed field is of 
the same type as the original one. Furthermore, the size parameter A' of the 
Nahm transformed field is given by: 

A' = -A . (63) 

Therefore, the Nahm transform provides a nonlinear relation for our pertur- 
bative expansion. A full analysis of this point is difficult and lengthy and 
will be left out from this paper, however it is instructive to look at the first 
few terms of this connection. 

In order to construct the Nahm transform one has to study the zero modes 
of the Weyl equation in the fundamental representation of the group: 

{V - 2mz) z) = , (64) 

where 2^ represent the coordinates of a point in the Nahm dual torus. The 
Weyl operator T> contains the self-dual potential Eq. (^ which can be ex- 
panded in powers of V AA. Although, for A 7^ 1 the configuration is not 
self-dual, it is still possible to define a Nahm transform (which will not be 
self-dual). Thus, we can expand xi^'-, -2) in the same way and equate to zero 
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all of the powers of the equation separately. From Eqs. (PTD,(P^) it is easy 
to see that the upper and lower components in colour space only mix for 
odd-even or even-odd powers of the expansion parameter. Thus, we might 
write: 

X{x;z)=\ - ■ '"^^'^ I . (65) 



In the previous formula, the explicit vector is in colour space, while the 
quantities 0±(x; z), (j)'^{x; z) are bi-spinors, which can be expanded in power 
series in AA. Eq. ( |64|) amounts for 0+ and (f)'_ to the equations: 



{Di - iS^ - 2Tiiz)(j)+ = iVXAW}^'__ (66) 
(D_i + iS^ - 2mz)(f)'^ = i^W^+ (67) 

and a similar equation holds for the remaining components. The symbol Di 
is defined in (^21). We see that in this way we get two independent solutions 
of (0) as predicted by the index theorem. The Nahm transformed SU{2) 
vector potential is then given by the formula: 

A'^iz) =ijd'x x'\x- z)^x\x; z) , (68) 

where the indices i,j G {1,2} label the two linearly independent and or- 
thonormal solutions. 

For the whole construction, the question of the boundary conditions sat- 
isfied by the spinors is crucial. Indeed, the naive periodicity requirement: 

TT X 

4>±{x + e^; z) = exp{±z-n^^-p} (f)±{x; z) (69) 

is inconsistent. How to remedy this situation is what is studied in Ref. In 
the case at hand the easiest way out is to impose the periodicity requirement 
only for the Xq and Xi direction, while requiring only double period conditions 
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on the other two. In short, this is just rephcating the torus in the X3 and X2 
directions. Consistently the integration in ( |68D has to be performed in this 
larger torus. 

To illustrate the procedure we will explicitly work out the lowest order 



term 0^ , which satisfies: 



(0) 

, VVliH..vli DCXljlOiiCO. 

(0) 



(Di - 2mz) = (70) 



This is a modification of the general equation studied in the previous chapter 
for q = \- Hence, following the same steps as before and imposing the new 
boundary conditions, we arrive at a unique solution (up to a multiplicative 
constant): 

ct)f{x-z)=eMT^iz^y,)}'^{y) ^ (71) 
where we have defined the auxiliary variable: 

~ ~l~ '^l^luTT'^u^u ) (72) 

and the function is the same one that appears in expression (p5D, with the 
replacement of /2,3 by 2/2,3 ^-nd ro_i by 2tq i. The constant K'^'^'' is fixed by the 
normalisation condition. In the same way one derives for 0?-' the expression: 

cl)^^\x;z)=tT2(Pf^*{x;-z) . (73) 

Now we might compute Aj} (x) by replacing Eq. (^Tj) in (^8|) . Now intro- 
ducing the complex variables: 

Vt, = ^{yf, + in^^y^) (74) 

we can express all derivatives with respect to 2^ in terms of derivatives with 
respect to Vq, Vi, Vq and f^. For example: 

^ ^(0) _ f^... O.J f ^ ^ \\ AO) 



^..yo-2>;3^--^||^r. (75) 
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has a very simple dependence on Vq and hence one has: 
d (0) _ J f TT in \ (0) 



The result of differentiating with respect to Vq is much more complicated, 
involving derivatives of Riemann's theta function. However, in the expression 
for A it is possible to integrate by parts and make the derivatives with respect 
to Vq act onto the complex conjugate of for which the complex conjugate 
of Eq. (|76D allows us to obtain a simple expression. We end up with: 

(A^^Yi^) = ^jd'x 0f ^ {niyo + Aihk^ |-^t;o + y (^o + ^^3)}) 

= 27r/o/3^3 • (77) 

In the right hand side of the first equality we have displayed the contribution 
of the two terms entering the right hand side of (ffSl) . One can compute the 
other components of the Nahm transformed field in the same fashion arriving 
at: 

Af{z) = -B'^{z) (78) 

where B'^{z) is given by the same expression as in section 2, but with 
the lengths of the torus replaced by those of the Nahm dual one i'fj, = 
The previous result implies that the Nahm dual of constant field strength 
configuration is a constant field strength configuration, even if they are not 



self-dual, generalising the result of Ref. 22 . 



5 Conclusions 

In the previous sections we have presented a systematic expansion which 
allows the construction of self-dual SU{2) Yang-Mills solutions with twist 
tensor ^03 = "^12 = 1 on the torus. The size of higher order corrections 
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depends on the lengths of the torus. For certain torus sizes the solution 
becomes equal to the the well-known constant field-strength ones. The mag- 
nitude of higher order corrections grows as we move away from these torus 
sizes. We have also compared the landscape of the solution, as obtained from 
our analytical expressions to leading non-trivial order in the expansion, to the 
numerical result obtained through standard techniques. The result is quite 
satisfactory at a qualitative and quantitative level. Finally, the expansion is 
used to obtain the Nahm transform of the self-dual configuration. Curiously, 
the perturbative construction of the Nahm transform has the same structure 
as the direct perturbative construction of the self-dual configuration on the 
Nahm-self-dual torus. We have not been able to equate these expansions 
order by order, but have shown this to be the case to lowest order. 

Let us now comment about the usefulness of our programme. Up to a 
proof, which we do not give, of convergence of our expansion, our method 
gives a direct proof of the existence of the solutions and for them to have the 
correct number of degrees of freedom. Even if calculating higher orders of the 
expansion turns out to be a difficult task, the expansion can be of theoretical 
interest for different reasons. For example, it might allow to investigate some 
general properties of the solutions. The case of the interplay with the Nahm 
transform is interesting, and should be pursued. Furthermore, one can study 
certain extreme limits of the torus sizes, which might allow to obtain exact 
solutions. 

Finally, we comment that our choice of SU (2) and of the aforementioned 
twist tensor has been dictated by simplicity. There is, however, no a priori 
essential difficulty in generahsing the construction given here to other SU{N) 
groups and different twist tensors. 
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Figure 1: Invariant densities $33 for the configurations A, B, C of Table |I|, 
from top to bottom, are shown. Plots in the left column display the ana- 
lytical perturbative result, and plots in the right column display the exact 
(numerical) result. 
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Figure 2: Invariant densities $-[1 for the configurations A, B, C of Table |I|, 
from top to bottom, are shown. Plots in the left column display the ana- 
lytical perturbative result, and plots in the right column display the exact 
(numerical) result. 
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Figure 3: Invariant densities $23 ^^e configurations A, B, C of Table |I|, 
from top to bottom, are shown. Plots in the left column display the ana- 
lytical perturbative result, and plots in the right column display the exact 
(numerical) result. 
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Figure 4: Comparison of exact and perturbative results for the electric field 
Euclidian time profiles $33'' (left column) and (right column) in the 
configurations A, B, C, from top to bottom, of Table |l]. The Xq coordinate 
has been rescaled with the lattice spacing a as defined in the text, and the 
center of the object has been set as origin of coordinates. The interpolating 
lines are intended to guide the eye. 
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